Heat conduction and long-range spatial correlation in ID models 
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Heat conduction in one- dimensional (ID) systems is studied based on an analytical S-matrix 
method, which is developed in the mesoscopic electronic transport theory and molecular dynamic 
(MD) simulations. It is found that heat conduction in these systems is related to spatial correlation 
of particle motions. Randomizations of scatterers is found to break the correlation, hence results 
in normal thermal conduction. Our MD simulations are in agreement with the theoretical expec- 
tations. The results are useful for an understanding of the relation between heat conduction and 
dynamic instablities or other random behavior in ID systems. 
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I. INTRODUCTION 

The study of heat conduction in one-dimensional (ID) 
systems is an interesting subject in the context of 
nonequilibrium statistical physics, which has been at- 
tracting much attention in recent years [1]. Most works 
aim at understanding the dynamic properties of heat 
transport in ID systems [2-10]. Although many mod- 
elling systems have been carefully studied and their ther- 
mal conductivities and temperature profiles calculated, 
the key dynamic properties of normal conduction is still 
unknown. Generally, normal heat conduction is specified 
by a diffuse-type motion, hence the investigation of ran- 
dom behavior of systems is necessary. It is well known 
that chaos can generate the needed random behavior in 
some systems. For example, thermal conductivity is char- 
acterized normal in the ding-a-dong model [2] and the 
Lorentz gas channels [6], due to their dynamic instabil- 
ity. However heat conduction is still considered abnor- 
mal in other systems, such as Fermi-Pasta-Ulam (FPU) 
model [3] in spite of their dynamic instability. Recently, 
Li et al. [9] studied three models with the same zero- 
Lyapunov exponents and found that heat conduction can 
be either normal or abnormal depending on details of the 
models. However they could not reach the conclusion 
that there is direct connection between dynamic insta- 
bility and normal heat conduction. From these studies, 
the most important and general conclusion is that the 
momentum conservation systems with non-zero presure 
have anomalous conductivity, whereas, anomalous con- 
ductivity does not always imply momentum conservation. 
Prosen and Campbell [11] were the first to prove this con- 
clusion. Recently, Narayan and Ramaswamy [12] criti- 
cized that the original proof is questionable, but there is 
no strong objection to deny the conclusion. The general 
theoretical result will explain the abnormal heat trans- 
port found in FPU model, and may clarify the recent 
confused results obtained by numerical works in some 
models, such as diatomic Toda model. Heat conduction 
in the Toda model was thought as normal by Jackson et 
al. [13] and by Garrido et al. [14], but was against the 
numerical results by Hatano [5], Dhar [8], Savin et al. [15] 



and Grassberger et al. [10]. This indicates that it is dif- 
ficult to get definite conclusions about the macroscopic 
transport properties of nonlinear systems by merely using 
short-time numerical simulations in small systems. 

In this paper, based on both theoretical and numerical 
methods, we try to explain the required dynamic char- 
acteristics that guarantee normal heat conduction in ID 
systems. By considering a ID model with N scatter- 
ers and noninteracting classic particles, we find that the 
classic model can be theoretically treated using S-matrix 
theory [16] developed in the mesoscopic electronic trans- 
port theory (METT). In METT, if the scattering is ran- 
dom, electronic transport will be incoherent, electronic 
phase correlation at the two ends of a ID chain will be 
absent, and eventually the Ohm's law is observed. Simi- 
larly, in our model, defining a 'phase' correlation, we find 
that normal heat conduction is characterized by the bro- 
ken of phase correlation, where the randomization of the 
scatterers makes a main contribution. This conclusion is 
supported by some principle considerations and previous 
related researches and can be applied to more general 
systems. Our result indicates that a sufficient random 
characteristic of ID systems is necessary to obtain nor- 
mal heat conduction. The random may be ensured by 
dynamic instablity, randomization of scatterers or other 
matters and could destroy long range correlations of par- 
ticle motion along the ID chains. Therefore, dynamic 
instability is not a necessary condition to guarantee nor- 
mal transport if other random effects exist. 



II. THEORY 

A. Model and Scattering Theory 

Consider a ID classical chain with N scatterers s^, 
{k — 1, • • •, A^) and non-interactive particles, which are 
elastically transmitted or reflected at these scatterers, we 
place this chain between two heat baths. Without loss 
of generality, the length of chain (L) is set as A^. If the 
transmission coefficient (t^) of each scatterer is a num- 
ber between and 1, it is a typical random walk pro- 
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cess, hence heat conduction is normal and independent 
of the kind of heat baths. However, here we suppose that 
scatterers periodicaly turn on or off in time, so parti- 
cles will completely transmit through a scatterer in some 
time ranges, but reflect at the scatterer in the other time 
ranges. So the transmission coefficients are function of 
time 4>, {tk{(l)))- It is a typically deterministic system. 
Here, we set the period to 1 and the average trans- 
mit coefficient to 1/2. We suppose that there arc different 
initial time shifts 5k for different scatterers, then we have 
tfe(^) = t{(l)-5k), and 



^^'^^ 10, l/2<</.<l. 



(1) 



The term (f) as well as 5k are thought as 'phases'. Another 
parameter of the fcth scatterer is its positions Xk- {5k} 
and {xk} will determine the properties of the model. Ac- 
tually, the model is very similar to the Lorentz channel 
model [6] or the Ehrenfest gas channel model [9]. 

We define the average of the particle current density 
J{x, V, t) on the model. 



M-l 



j{x,v,<j);n) = J{x,v,(j) + m + n), 



(2) 



m=0 



where n, m are integer numbers, M is a large integer 
number, and < 4> < 1. In steady state, j{x,v, 4>;n) is 
independent of n, and is noted as j{x,v,<j)). Due to the 
particle current conservation, if x is between two nearest 
neighbor scatterers, Xk and Xfe+i, J{x,v,t) is written as 
Jy{t — x/v). Hence we have, 

j{x, v,4>)='^ j^\v) exp {2m7ri {x/v - (/))}, (3) 

m 

if Xk < X < Xk+i. 

The heat current density Ju{x) and temperature profile 
T{x) can be writen as. 



J, 



and 
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'u{x) = ^j^ dvv^\j\;;\v)-jl^\~v)], (4) 



T{x) = 



(5) 



respectively, while Xk < x < Xk+i- 

Based on the properties of scatterers, we easily obtain 
the following scattering formula. 



P(k) 



(6) 



where £, f, and r''^^^ are nothing but S- matrix el- 

ements. Since this is a multi-mode scattering problem, 
each S-matrix element is an infinite dimensional matrix. 



For example, the element t[nri of fC^) means a coefficient 
which the mode m at the left transmit the scatterer k to 
mode n at the right, 

ttiiv) = t^""-"^ exp{-27ri(m - n){xk/v - 5k)}, (7) 

where t'^P^ = d4> t(0)e2P'^*'^, is the pth Fourier's ex- 
panded coefficient of transmission function t{4i). Simi- 
larly, other matrixes can be writen as, 

rWL{v) = ffe"""^ exp[-27ri(m -|- n)xklv\, 



Crf C^') = i^"""^ exp[-27ri(m - n){xkl{-v) - 5k)], (8) 
r'^^ {v) = ri'"-"^ exp[-27ri(m + n)xk/{-v)], 

where r^^'' = {5p^o — t^^^) exp(2p7ri 5k)- 
For whole N scatterers, we have. 



j(^\v) \ _ (f {v) R'{v)\ f J^Kv) 



R{v) f'{v)J \j^''\-v) 



(9) 



where f, R, T' and R! are whole transmitting and re- 
flecting matrix of all scatterers from left to right and from 
right to left, respectively, j^^-' and j*^^' correspond to the 
current densities at the left and right ends, respectively. 

Since the parameters of the system [xk and 5k ) deter- 
mine the S-matrix, we study heat conduction of models 
with different parameters. Without loss of generality, we 
set 



4 = c * i?, 

Xk = k-{).^ + d*{R-Q.b), 



(10) 
(11) 



where i? is a random number uniformly distributed be- 
tween and 1, c and d are the magnitude of disorder in 

scattering phases and positions, respectively. If both c 
and d equal to zero, it is a periodic scattering system, 
otherwise, it is a disordered system. For the latter, we 
can consider two kinds of different disordered systems: 
(1) the dynamic random system (DRS), where Xk or 5k 
are random in time; (2) the static random system(SRS), 
where random Xk or 5k are fixed rather than depend- 
ing of time in each realizations of system. We are in- 
terested in the average properties of many realizations. 
DRS and SRS correspond to the phonon scattering and 
impure scattering in electronic transport, respectively. 
Obviously, the disorder of Xk or 5k will induce some ran- 
dom phases in the S-matrix, it may affect the transport 
property of the system. 

The combining S-matrices of TV scatterers can be eas- 
ily written according to the well known rules generated 
in METT: it is written as a summation of many terms, 
where each term is depicted by a 'Feynman path' [16]. 
For example, the combining transmitting matrix of any 
two parts of scattering s^^^ and s'^' is. 



£(12)=f(2)[7_/.'(i)^(2)]-if(i)^ 

= f(2) [7 + ^'(1)^(2) +...]£(! 



(12) 
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Generally, wc have Tmn{N) = Ap, where each is 
a complex number contributed from path P which starts 
from mode n at the left end, to mode m at the right end. 
If tk ((/>) of each scatterer is independent of time </>, then 
there is only one mode, T = J2po ^^Poj where Pq are all 
the 'Feynman' paths with mode for all scatterers. We 
easily obtain, 



m) = i-^ + 1) 



N{l-t)+t' 



(13) 



While ^ oo, we have T{N) ~ 1/N, satisfying the 
Fourier's law, which is the expectation of the random 
walk process. 



B. Heat Baths 

Before going on to the study of the transport proper- 
ties of the model, wc analyze the effects of heat baths. 
Here, the temperatures of the left and right heat baths 
are noted as Ti and T2, respectively. Many kinds of 
heat baths can be selected. For example, we can choose 
complete-reflecting heat baths: as a particle hits a heat 
bath, it will be reflected with a velocity distribution 
Pt{v), where T is the temperature of the heat bath. 
Pt{v) can be of the Maxwellian distribution, Pt{v) = 
v/T exp(-?;V2T), or of a single velocity distribution, 
Pt{v) = 5{v — Vt), or other forms. In this case, we have 
completely reflecting conditions at the boundaries: 



(14) 
(15) 



where v is larger than 0, jln{v) / jm{~^) corresponds 
to the currents entered into the ID system from the 
left/right heat bath {x = 0/L), which is a different 
phase from j^{v) / jmi^"") which are used in eq. (9). 
•^i = '^'"^^rni-v) and j;, = dvf^{v). Actually, 
and are the Fourier expanded coefficients of the 
current densities J'((^) (at the left end, x = Q) and J''((^) 
(at the right end, x = L) oi ID chain, respectively. Prom 
eq. (9), we easily obtain, 



frnniv) Knniv) 
Rmniv) f^nip) 



(16) 



where there are some phase differences between the ma- 
trix elements and the original S-matrix elements in eq. 
(9), 

Tmniv) = exp(2rmTiL/v) T^n{v) 

TLniv) = T'^n exp(2n7riL/t;) (17) 

R'^^{v) = ex.p{2rmriL/ v) R'^n exp(2n7rii/t;) 

Then, we know and are not arbitrary, but satis- 
fying the following equation, 



^ Rmn ^ 1 



< T' 

mn 



>i 



' -It™ 



>2 
>2 -S„ 



(18) 



where < /„„ {i = 1, or 2 ), means the average 
value of fmn{v) under the velocity distribution function 
Ppiiv). The physical meaning of eq.(18) will be easily 
understood. In fact, we can treat the right boundary as 
{N + l)th scatterer ( complete reflecting), must sat- 
isfy the linear equation, (i?^'' — = 0, where R^^ is 
the total reflecting matrix of all N scatterers as well as 
the right heat bath. Prom Eq. (18), we have 



Nr 



R 



Since 



--<R>i + <f' >2 [/- < R' >2]~^ < r >i . (19) 



[I-R]- 



I + R + RR + - 



(20) 



we find that the result is nothing but the summation of 
'Feynman path'. Similarly, we also have the formula of 
i?'^ and equation (i?'^ — I)J^ = 0. If there is only one 
linear independent solution of the equation, we will get 
a unique steady state, (a free constant Jq is decided by 
the particle density and average temperature of system). 

However, the asymptotic N dependence of the heat 
current is only determined by the dependence of trans- 
mission coefficient |T|, which is independent of the details 
of heat baths. In this paper, we only consider a simpler 
case, replacing with the completely reflecting heat baths: 
particles uniformly enters into the system from two heat 
baths at time (</>) with the single velocity distribution 
6{v — Vt), therefore the heat current is simply propor- 
tional to the transmission coefficient Too. When calcu- 
lating T{x), we choose another uniform incoming current 
from the right heat bath at the same time, making the 
total particle current zero. 



C. Phase Correlation 

It is well established that, if the scattering is random, 
the contributions from different Feynman paths are in- 
coherent in electronic transport, therefore, we can find 
ohm's law and normal electric conductivity. But if there 
are some long range correlations between scatterers, the 
electronic phase-relaxation length may be larger than the 
length of the system, the transport will be coherent and 
we can find anomalous electric conductivity. So the ex- 
istence of the phase correlation between electrons at the 
two ends can be used to judge whether the transport be 
coherent or not, hence whether the conductivity be ab- 
normal or normal. Comparing the results, we expect a 
similar relationship between heat conduction and a long 
range correlation of particles in ID systems. We also 
expect the random charactaristic of scatterers to be re- 
sponsible for the correlations. 

For any incoming particle current from the left heat 
bath, J\.{4>), we have the transmit current density at 
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the right boundary J^{tjj) = J T{tl},(j))j''_^{(j))d(j), where 
T('i/;, 0) is the transmit function, 

T(V', (t>)=^ cxp{-2m'Kiip) Tmn exp(2n7ri0) (21) 

mn 

Similarly, the reflecting function is -R(V') </>)• Since for any 

(t>, J[T{i), 4>) + R{^, 4>)]dtP = 1, so Ron = Son - Un- 

The probability that we observe a current with phase 
tp at the right end and a current with (j) at the left end 
is W2{'ip,L;(j),0) = T{tp,(j))j\^{(j)). We define a normal- 
ized two-point current distribution function as /2(V') 4') = 
W2{i^, N; (j), 0)/A, where A is a normalized constant, 

A = j W2{'>p,L-(t>,Q)d(j)di) = j r+{i))di). (22) 

The left and right normalized distributions are 



and 



= J /2(V', 



(23) 



(24) 



respectively. 

Based on these distributions, we define a motion cor- 
relation of particles at two ends of the ID chain as. 



D 



W)2 - 



aidr 



(25) 



where (• • •) ■ (i = 2, ^ or r) means the average value under 
the distribution functions f2{il>,^), fr{'4') aiid fi{^), re- 



spectively, ai/r = \J {{4>^)i/r ~ {4>)i/r) ^hc width of the 
current distribution at the left/right end. Hence, D = Q 
implies the broken of the correlation and T{ip,(j)) = 
F{tJj)G{4') (or Tmn = FmGn)- Here, D corresponds to 
the electronic phase correlation at both ends of ID elec- 
tronic transport system, and the particle velocity play 
the role of wavevector of electron. 



D. Theoretical Results 

Based on the rules of the 'Feynman' paths on scattering 
theory, the transmission matrix T^n can be writen as, 



Tmn — tmptjn -^pqi^)' 



(26) 



where f*^^^ and t^^^ are the transmission matrix of the 
scatterer N and the scatterer 1, respectively. Bpq(N) is 
the sum of all 'Feynman' paths which start from mode 
q of the scatterer 1, end to mode p of the scatterer N. 
Expanding Bpq{N) in large N, the lead term is noted 
as apq/N^ , therefore we have the first result that the 



asymptotic A'' dependence of T^n is independent on m 

and n, 



Tmn = T{N)h 



(27) 



where T{N) ^ 1/N^, is a simple notation of Tqq and 
hmn is independent of N. If the correlation is absent, 

Tmn = T{N)fm9n. 

For DRS, due to the random time-dependent posi- 
tions/phase of the scatterers, the contribution from dif- 
ferent Feynman paths are not coherent with each other, 
so they can be first averaged in time, then be summed 
up. This indicates that only 0-mode contribute to the to- 
tal transmission coefficient, so the heat conduction shall 
be normal. We also easily know that the correlation is 
absent in DRS due to the incoherent Feynman contribu- 
tion. 



III. NUMERICAL SIMULATION 

In this section, by using a uniform inputing current 

with single velocity vi = \/T\ from the left end, we nu- 
merically simulate the transmit coefficient T(N), the dis- 
tribution function frijp), fii(t>) and the temperature pro- 
file T{x) in different systems. First, we find that both 
fr{ip) and fi(ijj) are independent of A'' for all models in- 
cluding the periodic scattering system, DRS and SRS 
with the same disordered magnitude (c and d). The 
results verify very well our theoretical expectations in 
eq. (27) and indicate that hmn is only dependent on the 
disordered magnitude. As c (or d) increases, /r(V') and 
fi{^) are flatter, but they are not uniform even though 
the disordered magnitudes arrive at their maximum value 
(1). For example, the results in DRS are shown in Fig. 
1, where J^{ip) is the unnormalized current distribution 
function at the right end. Then, we show the temper- 
ature profiles T{X/N) of DRS in Fig. 2(a), which are 
found to being independent on N but dependent on the 
disorder magnitude c (or d). It indicates that dT/dX 
is proportional to 1/N. The heat current J„ ~ 
is shown in Fig. 2(b), and fitted 7 w 0.9989 ± 0.0045. 
Therefore the heat conduction is normal in DRS with 
disordered scattering phases, which verifies our theoreti- 
cal expectation. In this case, our obtained spatial corre- 
lation D is very small but with the same order statistic 
error. For a better correlation estimation, we used the 
non- uniform input currents J\.{4') with different distri- 
butions and found the right-end distribution functions 
f^iip) independent of the selection of J\-{4>) (not shown). 
This indicates that the transmission function T{ip,(j)) can 
be writen as F{'tp)G{(p), hence D = 0. For DRS with 
disordered scattering positions, the obtained results are 
similar as that of DRS with disordered phases. 

We also numerically simulate the heat conduction of 
SRS. Since the position disordered system is similar to 
the phase disordered system, we only show the results 
of the latter. The temperature profile T{x) and the N 
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dependence of J„ in SRS are similar as that in DRS, but 
7 ft! 1.19 in SRS, slightly larger than 1 in DRS (seeing 
Fig. 3) and there is a slight difference on T{x/N) of 
systems with different A'^. The differences may be due 
to the limiting chain length N in our simulation. For 
different disordered magnitude c, the value of 7 is found 
to be indistinguishable. For comparison, we simulate the 
heat current J„ in a periodic scattering system, the result 
is also not very clear. It may be necessary to simulate 
them in longer chains. In Fig. 4, we show the distribu- 
tion of D for different disordered realizations in SRS. The 
average of D is about zero, in agreement with our expec- 
tation which the scattering disorders induce the absence 
of the correlation. Other interesting result is shown in 
Fig. 4(b): the fluctuation of D decreases as N increases, 
{D'^) = 7V-°-^2 indicates that D ^ in any single 
disordered realization of SRS while N ^ 00. Contrarily, 
in periodic system, the correlation D is not zero. 



variant of ding-a-dong model, which its even-numbered 
particles are oscillatedly coupled each other, rather than 
coupled with its individual lattice site in the initial 
model [2] . Due to the existence of some long- wave modes 
in the model, the motion correlation is obviously present, 
we expect the thermal conductivity should be abnormal. 
On general ID systems, important questions are how to 
define suitable correlations and judge directly the exis- 
tence of the correlations from the characteristic of sys- 
tems. For example, for momentum conservation systems, 
which thermal conductivities are anomalous if presure is 
not zero [11], it is worthy for connecting the conditions 
with existence of spatial correlations. 

X. Z. was supported by the Grants-in-Aid for 
Scientific Research of JSPS, his current address is 
zhou@pe.titech.ac.jp. X. Z. thanks the help of Mr. C.-Q. 
Li. 



IV. DISCUSSIONS AND CONCLUSION 



Wc have defined a spatical correlation along the ID 
chains to describe scattering (i. e. transport) charac- 
teristic in the systems and theoretically expected that 
randomization of scattering detcrminced the existence of 
the long range spatial correlation, therefore properties 
of heat conduction. For DRS, the expectation was con- 
firmed by both S-matrix theory and molecular dynamic 
simulation. For SRS and periodic system, according to 
the obtained correlations, normal and abnormal heat con- 
ductions were expected, respectively. But our numerical 
simulations in some small systems could not completely 
clear the expectation, more simulations might be neces- 
sary. A noted fact is that our model is very similar to 
the Ehrenfest gas channel [9]. In the latter, the channel 
is quisal-lD systems with a small transverse coordinate 
(the height of the channel). Actually, the height corre- 
sponds to the 'phase' in our models, a similar scattering 
theory can be derived. For the Ehrenfest channel, the 
results of Li et al. [9] , which heat conduction is normal if 
the heights or positions of the scatterers are random and 
it is anormalous if their heights or positions are periodic, 
can be easily understood from our results. In these cases, 
the required incoherent characteristic of motion to get 
normal heat conduction can be ensured by the random 
'phases' and positions of scatterers, not needing additive 
dynamic chaos (non-zero Lyapunov Exponent). There- 
fore, the conditions with normal heat conduction in ID 
systems are whether or not there are sufficient random 
to dcstory completely the correlation of motions along 
chains. The random may come from dynamic instablity, 
random scatterers or another sources, but no any partic- 
ular random, such as chaos, is necessary. The judgement 
of the incoherent motions (i. e. normal transport) is 
whether or not the long range spatial correlations exist. 

We also can consider another models, for example, an 
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X. Zhou et al. Fig. 1 X- Zhou et al, Fig.2 DRS (d=0) 




FIG. 1. The particle current distribution J^{(j)) at the right 
end in DRS with random phases: (a) J^{4>) of systems with 
different A'^ are identified each other, where c — 0.8. (b) the 
distribution {(j>) of DRS with different c. c = means a 
periodic system. 



FIG. 2. The results of DRSs with random phase: (a) Tem- 
perature profile, where x — X/N . For every c, we show 
four temperature profiles, the corresponding chain length are 
32, 64, 128 and 256, respectively. They are almost identified 
each other, (b) Heat current J versus A*'. 
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X. Zhou et al. Fig. 3, SRS cl=0 X. Zhou et al. Fig. 4 




FIG. 3. The results of SRSs with random phase and peri- 
odic systems: (a) Temperature profile of four chains, where FIG. 4. The correlation of SRSs with random phase: 
c = 0.8. (b) Heat current J versus A'', from top to bottom, c (a) The correlation distribution of different samples, where 
are 0.0 (periodic system), 0.2, 0.4, 0.6 and 0.8, respectively. A'' = 32 and c = 0.2. (b) The fluctuation of correlation (D^) 

versus A*'. We show the results of three SRSs with different c 
(c = 0.2, 0.5 and 0.8), they are almost indistinguishable. The 
solid line is the best fit one. 
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